Generalized boundary conditions for the circuit theory of spin transport 
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The circuit theory of mesoscopic transport provides a unified framework of spin-dependent or 
superconductivity-related phenomena. We extend this theory to hybrid systems of normal metals, 
ferromagnets and superconductors. Our main results is an expression of the current through an 
arbitrary contact between two general isotropic nodes. In certain cases (weak ferromagnet and 
magnetic tunnel junction) we derive transparent and simple results for transport properties. 



I. INTRODUCTION 

Spin-transport in hybrid ferromagnet (F)-normal 
metal (N) systems has been object of extensive inves- 
tigation since the discovery of the giant magnetoresis- 
tance (GMR) effect.EJ Magneto-electronic multiterminal 
systems lead to novel applications, e. g. magneto record- 
ing heads or magnetic-field based sensor devices. Fu- 
ture prospectives of non- volatile electronics motivate also 
many fundamental and applied research. Transistor- 
like effects were found in a ferromagnet-normal metal 
three-terminal device, that depend on the relative oiL 
entation of the magnetization of the ferromagnets Btil 
The dependence of transport properties with the rela- 
tive angle between the magnetization directions of the 
ferromagnets, was addressed experimentallyu as well as 
theoretically^ in the past. These days non-collinear spin 
transport attracts an increased attention due to the re- 
cent interest on the spin-current induced magnetization 
torques.ETtj A ferromagnetic single-electron transistor in 
a three terminal configuration has been realized and stud- 
ied theoretically.EE3 

If a normal metal (N) is attached to a superconductor 
(S), quasiparticles of different spins are coupled via Ait. 
dreev reflection on the normal side of the NS interfaced 
A strongly modified density of states (DOS) in the nor- 
mal metal caused,-by induced superconducting corre- 
lations was foundB This is so-called superconducting 
proximity effect. In contrast with normal metals, in a 
ferromagnet (F) the presence of a strong exchange field 
leads to big differences between the two spin bands. The 
question of the coexistence of ferromagnetism and super- 
conductivity has been extensively investigated in the last 
decades. The effect of spin splitting on the superconduct- 
ing -proximity effect was already investigated long time 
agctj. New experimental developments on proximity ef- 
fect in hallistic ferromagnetic layers have_been recently 
reportecO and theoretically confirmed.lij Results ob- 
tained in ferromagnet-superconductor nanocontacts have 
been explained in terms of bias dependent transparency 
of the FS interface .EjO Diffusive heterostructures of fer- 
romagnets and superconductors showed unexpected large 



values of the conductance of the ferromagnetic part. I 




Strong mutual influence between the superconductor and 
the ferromagnet and long range proximity-effects have 
been proposed to explain these rcsults.c3l2j Effects re- 
lated to the interplay between spin accumulation and 
Andreev reflection has been also investigated.^ On the 
other hand, the prediction of an exotic superconduct- 
ing state formed in a superconductor by the presence 
of an exchange field h, was independently reported b« 
Larkin-Ovchinnikov and Fulde-Ferrel, some years ago.E3 
Experimental confirmation of these effect in bulk super- 
conductors is still needed. Many investigations have fo- 
cused on the. study of thermodynamic properties of FS 
multilayers. 

Non-equilibrium Keldysh Green's functions in 
the quasiclassical approximation, have been exten- 
sively used in th£._uast to describe non-equilibrium 
superconductivity.EElE!] This quasiclassical theory of su- 
perconductivity is based on semiclassical transport equa- 
tions for quasiparticles. Generally to solve these trans- 
port equations is technically difficult, however the final 
results are relative simple and clear— .The so-called "cir- 
cuit theory of Andreev renection"E3, was conceived as 
a generalization of the Kirchhoff's theory of electronic 
circuits, to simplify these equations into a handful of 
accessible rules. The circuit theory of mesoscopic trans- 
port was recently extended tGMjlescribe transport in non- 
collinear magnetic structures^]. Interesting phenomena 
like spin precession effects on the induced spin accumu- 
lation have been foundo. An important concept in the 
theory of non-collinear j-spin transport is the so-called 
mixing conductance G^c3 which is closely related to the 
spin-current induced magnetization torquepp and the 
Gilbert damping in thin ferromagnetic films .Ej Moreover 
a three terminal spin-transistor device has been . pro- 
posed in the framework of the circuit theory.oEj The 
main advantage of the circuit theory description is that 
it provides a simple approach based on spin-charge cur- 
rent conservation to calculate the transport properties 
of (multi-terminal) mesoscopic hybrid systems. This 
is achieved by mapping the concrete geometry onto a 
topologically equivalent circuit, represented by finite el- 
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ements. We note that this very same description allows 
to calculate a v Vifltjf of transport properties, such a& 
current statisticsE3'E2l, weak localization correctionsEH 
or transmission eigenvaluesc3. Consequently, it is worth- 
while to extend this general framework to the all possible 
combinations of heterostructures. 

So far, the developed circuit theory is suitable to de- 
scribe electron and spin transport in multiterminal hy- 
brid ferromagnet-normal metal FN or superconductor- 
normal metal SN systems. Obviously, an extension of the 
circuit theory to hybrid systems combining ferromagnets, 
normal metals and superconductors is necessary. This is 
done in the present paper. We derive a general expres- 
sion for an arbitrary contact, which, however, turns out 
to be a little bit unhandy. To obtain manageable expres- 
sions we present also approximate results for two special 
cases. Below we list all our results: 

• The general matrix current Eq.(|l4|) through an 
arbitrary spin-dependent connector. It requires 
the knowledge of the full scattering matrix (or 
transmission matrix) as consequence of the non- 
separable Spin- and Nambu structures. This in- 
hibits the transformation into the normal eigen- 
modes of the scattering problem. Eq. ( |l4| ) is there- 
fore mostly of numerical interest. 



• The matrix current Eq.(19) for a weakly spin- 
dependent contact. Here an expansion in terms of 
normal eigenmodes is possible. Eq. (19) is a spin- 
dependent correction to the matrix current from 
Ref. (||). 

• If the spin-dependent contact is a tu nnel barrier, 
the tunneling matrix current Eq. ( |21a| ) takes a par- 
ticular simple and transparent expression. The 
properties of the contact can be expressed in terms 
of the spin conductancesG^Q) and the (complex) 
mixing conductance G^Ej. 

The paper is organized as follows. In section || we illus- 
trate some basic aspects of the circuit theory. In Section 
[II we present the microscopic description of a contact re- 
gion or connector in the circuit theory. In section IV we 



calculate the matrix current through a general connec- 
tor between two metallic regions (nodes) of the circuit. 
The nodes can be of N-, F- or S-type and the contact 
is assumed to have a arbitrary magnetic structure (e. g. 
magnetic tunnels junctions, magnetic interfaces). This 
expression constitutes the generalized matrix current for 
the circuit theory. In section ^ we make use of a pertur- 
bation expansion in the spin structure to obtain simpli- 
fied expressions for two cases. Details of the pertubation 
expansion are presented in the Appendix. In Section VI 
we present our conclusions. 



II. CIRCUIT THEORY 

In the circuit theory the system is split up into reser- 
voirs (voltage sources), connectors (contacts, interfaces) 
and nodes (low resistance islands/ wires) in analogy to 
classical electric circuits. Both reservoirs and nodes are 
characterized by l6N c h x 16iV c /j Green's functions G, 
which are matrices in Keldysh<g)Nambu<g>SpinCg)Channels 
space, where N c h are the number of propagating modes. 
These Green's functions play the role of generalized po- 
tentials of the circuit theory. The Green's functions in 
the reservoirs and nodes are assumed to be isotropic in 
momentum space. This requires that sufficient elastic 
scattering is present both in reservoirs and nodes, due to 
the presence of random scatterers and irregularities in the 
shape. This justifies the use of the diffusion approxima- 
tion to describe transport. This assumption is reasonable 
since hybrid (multi-terminal) devices are quite dirty sys- 
tems. In general for the stationary case, G depends on 
space coordinates and energy 



6(^,6) = 

being 
G(ry,t-t') 



dt ',t - t') exp{i-(t - t')} (la) 



G R {r,r*,t-t') G K (r,r*,t-t') 







G A (ry,t-t>) 



(lb) 



G R (r,r*,t-t') = -i6(t-t')({i>(r),¥(r")}) (lc) 
G A (f,r<,t-t')=i6(t' -t)({$(r),¥(r*)}\ (Id) 



G K (r,r',t-t') = -i{ 



(le) 



where [.., ..] and {.., ..} denotes commutator and anticom- 
mutator respectively. Semiclassical and diffusive approx- 
imations allow to obtain Green's function, which depends 
only in one spatial coordinate G(r, e)o On the nodes 
the Green's functions are also assumed to be spatially 
homogeneous, depending only on energy G(s). This re- 
quires that the resistance of the node is much smaller 
than the contacts resistances connecting different nodes, 
which implies also that the current through the device is 
controlled by the contacts resistances. Regions with spa- 
tially dependent Green's functions (e.g. diffusive wires) 
are modelled by an appropriate discretization. For ex- 
ample, a quasi one-dimensional diffusive wire can be rep- 
resented by a series of tunnel junctions and nodes. In- 
ternal dynamics along the wire (finite energy transport, 
spin- flip, etc.) is inclu< 



in this description as a leakage 

current from the node.E3 

In analogy to Ohm's law in classical electric circuits 
(/ = V/R), in the present circuit theory it is essential to 
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obtain the "spin-charge" matrix current that flows be- 
tween two nodes (or between a reservoir and a node) 
through a contact/connector. This matrix current de- 
pends on the properties of the connector (analog to the 
resistance R) and on the Green's functions at both side 
of the connector (analog to the voltage drop through the 
connector V ). In general the matrix current I{z,e) is 
defined in Keldysh®Nambu®Spin£g)Channels space as 

-^a,,,,,,),,., (2) 

where ? = (fi,z), z being along the transport direction 
and p being perpendicular to the transport direction. 
The electric current I e is defined in terms of the matrix 
current I{z 1 e) as 

1 t°° Ac 

where a, f are Pauli matrices in spin and Nambu space. 
G(f, r*;e) can— be expanded in transverse modes in the 
following wayci 

G(r,r<;e)= E 6 nw(^')x (3) 

nm,(T,(r' 

exp(iak s n z ~ ia'kiz') X s n (p)xi*(p')- 

Here n (m) label the propagating modes, a, a 1 = ±1 are 
the direction of propagation, s, s/ = {t, 1} are spin in- 
dices, fc* is the longitudinal momentum and Xn(p) is the 
transverse wave function. By using this representation, 
the current can be written like 

r< s '(z,e)^ e 2 £ (av s n + a'vi)G s n ^ la ,(z,s) (4) 

n ; a.a' 

x j dp x s n(p, z)x s m *(p, z )- 

Here v'Jfy = h k'Jfy/ m is the velocity in mode n(m) 
with spin s(s'). The transverse wave functions XniP> z ) 
are eigenfunctions of the Spin operator S = h <j z /2 and 
are normalized in a way that the carry unity flux current. 
In this case the current reduces to 

I s > s (z,e) = 2 i e 2 J2 * < e) (5) 

n ; a 

Note that the current I s ' s (z, s) in Eq.(|J) is a matrix diag- 
onal in spin space. In the case of a ferromagnetic material 
Eq. ([5]) is valid if the magnetization of the ferromagnetic 
material is parallel to the spin quantization axis. In the 
case of non-collinear transport where there are two dif- 
ferent magnetizations defined in the system M\ and M2 , 
Eq. (|^) has to be properly transformed into the basis of 
the spin quantization axis 

p,s' = u fs,sjj-l 

U being the spin rotation matrix that transforms M 1 ( 2 ) 
into the spin quantization axis. Note that in this case the 
transformed matrix I s ' s can have general spin structure. 



III. ARBITRARY CONNECTOR 



The contacts or connectors between the reservoirs 
and nodes are described in terms of the scattering for- 
malism by a transfer matrix M. The goal is to ex- 
press the matrix current given by Eq.(^|) in terms of 
isotropic quasiclassical Green's functions Gi(2) & t both 
sides of the contact and transmission/reflection coeffi- 
cients that characterize the scattering in the contact 
region. M is in general a 16iV c ^ x 16N c h matrix in 
Keldysh<g>Nambu<g>Spin<g>Channels space. Particularly 
M is proportional to the unit matrix in Keldysh space 
and diagonal in Nambu space. This is denoted by (~). 

The Green's functions G = G^*,(z,z') is not 
a continuous function at z = z/OO The values of 
G^lJ ' ,(z, z') for z > z' and z < zi are matched at 
z = zi by using the following representation (for details 
see Ref . g|||) 

-S, S t f\ 

n<T,m(r' ( Z ,^)= == (6) 

Y u n Urn 

| Q-Saa' $nm sjgn(z - Z 1 ) - 

V n 

where the matrix g{z,z') = g^ mrT ,(z, z') is continuous 

at z = zi . Note that the matrix g^ a s m(T t has non-trivial 
structure in channels space. We call this function ballis- 
tic Green's function. At the contact region (z — 0) the 
transfer matrix M = M^' a s , connects g(z, z') at both 
sides (left z, zi = _ and right z, zi = + ) 

9(2) = M g {1) M f (7) 

where g {2 ) = K z = 0+ > z ' = 0+ ) and 9(1) = d( z = 
_ ,z' = _ ). Note that in Eq.(f?]) there is implicit a 
sumation over Channel and Spin indices. 

An important assumption of this theory is the 
isotropization assumption. The ballistic Green's func- 
tion g defined at both sides of the contact region becomes 
isotropic in momentum space, when departing from the 
contact region. Such isotropization happens in a region 
of the order of the elastic mean free path Z m f p but much 
smaller than the spin diffusion length l a { = y/D t s { and 
the coherence length of a superconductor £a = y/D/A. 
In the isotropization zone the dominant contribution to 
the self-energy is the elastic scattering, which implies suf- 
ficient disorder (or chaotic scattering) in this region. The 

v 1 ( 2 ) v v 

self-energy is then = -i G 1{2 )/2 T imp , G 1(2) being 

the isotropic quasiclassical Green's function for the left 
(1) and for the right (2) side of the contact region, which 
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is proportional to the unit matrix in Channels space and 
i~imp is the impurity scattering time. In order to assure 



that the G^ ma ,(z, z') does not diverge at z(z') — > 0, we 
have to impose the following conditions!! 



(£* + Gi) (£*-&) =0 
(E'+Si) (E z -G x )=0 
(E 2 -G 2 ) (E 2 +.g 2 ) =0 
(£ 2 - sa) (E« + G 2 ) = 0, 



(8) 



(9) 



(10) 



(11) 



E 2 =cr Saa' being the z-Pauli matrix in the direction of 
propagation "sub-space" {ere'} of the Channels space. 
An important consequence of this approach is that af- 
ter the isotropization zone, the ballistic Green's function 

9nama' ( z > z ) ec l ua l s the isotropic quasiclassical Green's 
function Gmv 



IV. GENERALIZED MATRIX CURRENT 

In previous work,0 M and Gu2) commute because no 
spin structure was considered. Now due to their general 
spin structure, M and G do not commute. We multiply 
Eq. (g) and Eq. © by M from the left and by from 
the right. Adding both resulting equations, using Eq.(0) 
and using the normalization condition G 2 = 1,lj we get 
the following expression for gi 

gi = (1 + GiMt G 2 M)^ 1 [2Gi + (I - GiM* G 2 M) E z 

(12) 



Note that in Eq.(12) all the structure in channels is con- 
tained in M. From Eq.(Q), the matrix current / is be 
expressed in terms of §1(2) aso 



1(e) = 2 ie 2 J2°Vn G s na % a ,(e) = e 2 Tr„, CT [t z ~ 9l ] 

n ; a 



(13) 



Note that the Green's function G^J n<7 , is assumed to be 
spatially homogeneous, depending only on energy G(e). 
In this case that current does also not depend on posi- 
tion. By using the cyclic property of the trace, we find 
finally for the matrix current /(e) 

1(e) = e 2 Tr„, CT {(l + G X M^ G 2 M)^ (l - G±M^ G 2 M) 
+2 (l + GiAf t G 2 M) 1 t z Gi}. 

Eq. ( |i"4| ) is the most general expression for the current / in 
terms of isotropic quasiclassical Green's functions Gi(2) 



at both sides of the contact and transmission/reflection 
coefficients of the contact region. Once such expression of 
1(e) is obtained, we can apply the generalized Kirchhoff's 
rules, imposing that the sum of all matrix currents into a 
node is zero. This completely determines all properties of 
our circuit. Eq.(p^) therefore completes our task to find 
the generalized boundary condition for the circuit theory. 
However, in this form a concrete implementation requires 
the knowledge of the full transfer matrix (or equivalently 
of the scattering matrix). Usually this information is 
not available for realistic interfaces and one tries to re- 
duce Eq. ( |l4| ) to simple expressions by some reasonable 
assumptions. In the next section we will do this for two 
special cases. Note, that for a spin-independent interface, 
Eq. ( |l4| ) can be expressed by the transmission eigenvalues 
only, which is a formidable simplification. This transfor- 
mation is not possible anymore for the spin-dependent 
contact. 



V. TWO SPECIAL CASES 

We want to obtain more transparent and clear analyt- 
ical expressions of Eq.(14). The price to pay for that is 
loss of generality, since we have to make certain assump- 
tions about the contact. The main difficulty lies in the 
inversion of the matrix 1 + G\M^GiM in Channel-Spin 
space. Let us assume that the transfer matrix M can be 
split in the following form 



M = M a + SM = M (i + 6X) . 



(15) 



where Mo is a transfer matrix with structure in channel 
space only and proportional to the unit matrix in spin 
space, while 5M = Mq5X includes non-trivial structure 
in spin space. Note that the matrices Mo, Gi( 2 ) commute 
with each other, whereas in general SX does not commute 
with Mq and Gi( 2 ). Assuming that SX <C 1, we can per- 
form a perturbation expansion of (I + G\M^ G 2 M) 
in the parameter SX (see Appendix). The matrix M = 

(Ad) s,s , can now be diagonalized in the basis of eigen- 

V / ncr,mcr 

modes N(N) : (M) s ' s ' , -> (\1)~ ~ N ,. 

To 0th order M = Mq commutes with Gu2y In this 
case Eq. ( |l4| ) reduces to 



f°\e) = e 2 Tr N , 



1 



1 + Q GiG 2 



(2 t z Gi + l-Q GiG 2 ) 
(16) 



The hermitianj-cnnjugate matrix Qq (Q q = Qo), reads in 
(iCjhannel spacec 2 ] 



Qo 



A B 
St A 



(17) 



being A real and B complex N c h x N c h matrices. 
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The eigenvalues of Qq, appear in inverse pairs 
{qntIn 1 ), and are related with_.the transmission coef- 
ficients TV in the following wayc3 

-l 



qN + q N 



2-T 



N 



N 



being 



\B\ 



A 2 



1. 



By performing the trace over directions of the mode in- 
dices cr(cr'), the expression for the Oth order current re- 
duces to 



[G 2 , G\ 



N 



A + T N ({Gi,G 2 }-2) 



(18) 



Eq. ( |18[ ) is the expression for the matrix current obtained 
in Ref.( |3^) for a spin-independent contact. 



A. CASE I: WEAK FERROMAGNETIC 
CONTACT. 

Now we concentrate in the first order term I^(e), 
given by Eq.( A4c[ ) of the Appendix, SX = A being 



anti-symmetric with respect to time-reversal transforma- 
tion. In particular (SX) . is anti-symmetric with respect 
to transformation over spin space and symmetric with 
respect to transformation over directional space. In this 
case (SX) . describes the case of a weak ferromagnetic 

contact (SX^ A ~ <t T3Y The unity vector M is in 

the direction of the magnetization, and <r, f are Pauli 
matrices in spin and Nambu space, respectively. On the 
other hand, the (SX) g may describe spin-flip processes 
due to spin-orbit interaction at the contact. At first or- 
der in SX the contribution given by (SXj s vanishes. In 
general, to threat spin-flip at the contact, we need to g o 
to higher orders in SX (SX 2 , SX 3 , ..) . From Eq.(|A3oj), 
Eq.(A31) and Eq.( |A32| ) (see Appendix), Jj can be writ- 
ten as 

2 



N 



T N ({G u G 2 }-2) 

[{tlSh + 8tlt 2 ,G 2 \+C,G 1 
2 

X 4 + T N ({Gi,G 2 }-2) 



(19) 



where 



C 



(4 5*2 



8t\t 2 )(- 



2T, 



N 



T, 



N 



(4 



St 2 + r-JSri 



(20) 

Eq.(|l9|) is the spin-dependent correction to the matrix 
current given in Eq. (|l8|) . This result constitutes an im- 
portant step in the applicability of the circuit theory for 
F|N|S systems, because it provides of a simple prescrip- 
tion to describe magnetically active interfaces/contacts. 



B. CASE II: TUNNELING BARRIERS 



For the case when the contact arc tunneling barriers 
(Tjy <C 1), is possible to neglect the term Jj yjT Gi, G 2 } — 
2) in the denominators in Eq.(18) and Eq. (|l9[) . Keeping 
terms of order Tn, t\5t 2 /T^, t[St 2 and r^Sri, we can 
write the total matrix current 1(e) in a very transparent 
way like 



1(e) G T * -I G M r 



{M£f 3 ,G 2 },Gi] (21a) 



+i- 



,G d 



M a fs, Gi 



where 



- N 



Gt — Gq 2J Ti 

N 

Gmr = Gq ^ STn 

N 

iG^/2 = G Q ^(4^2 + r\Sr x 

N 

'T N ~2 



(21b) 
(21c) 
(21d) 



2T ; 



N 



t\8t 2 -St\ t 2 \) 



Gq = e 2 /2nh being the conductance quantum and where 
the explicit form of SX = SX M a 73 is used. For a weak 
ferro magne ti c conta ct, the elements St (Sr) introduced in 
Eqs. (Al8a)-(Al8d) in the Appendix can be expressed in 
terms of spin dependent amplitudes like 



(22) 



■ t TT 







' t + St 








t ll 







t-St 


" r TT 







r + Sr 















r — Sr 



(23) 



From these it easy to see that 
Gt = ^^GqT n = Gq 



1 n < 1 n _ fr "i ^4) 



N 



N 



Gmr — 2_, Gq ST n — G c 



rp\ rpi s~i\ S~tl 

N N T 



iV 



AT 



(25) 



Now we see that Gt is the usual conductance of the con- 
tact G T = (G T + G T ) /2, and G M r = (g t - G T ) /2 
accounts for the different conductances for different spin 
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directions, leading to a spin polarized current through 
the junction. 

On the other hand, for the case Tjv = (the case of 
a ferromagnetic insulator contact), Eq.(21a) is not zero 
and reduces to 



2 



M a f 3 ,Gi 



(26) 



where G^ depends only on reflection amplitudes at the 
normal metal side 



iG^/2 = G Q Y J 



r TT III 



t 



JV 



(27) 



A' 



The physical meaning of this term in this particular case, 
can be now understood as follows: electrons with dif- 
ferent spin directions pick up different phases when re- 
flecting at the ferromagnetic insulator. During this pro- 
cess, ferromagnetic correlations are induced in the nor- 
mal metal node. In this particular case, the coefficient 
G$ is related to the mix ing conductance introduced 
via G$ — ImG^. Eq.(21a) has been recently used as 
boundary condition current for the circuit theory applied 
in a MI|N|S structure, being MI a magnetic insulator 
and under conditions of superconducting proximity ef- 
fect. In Ref.( we shown that in this case, the effect of 
the conductance G^ can be seen as an induced magnetic 
field which give rise to spin-splitting of the induced "BCS 
like" density of states h = G < i > S/2Gq, being S the average 
level spacing in the normal node. In particular, we also 
show that for a system composed by two coupled MI|N|F 
trilayer structures the absolute spin-naive effect can be 
achieved for a finite range of voltages. □ 



VI. CONCLUSIONS 



which describe the adjacent reservoirs/nodes and trans- 
mission/reflection coefficients that characterize the scat- 
tering in the contact region. This expressions should be 
numerically implemented in order to solve any general 
arrangement of reservoirs, contact and nodes. 

Moreover, we have perform a perturbation expansion 
in the spin asymmetry of the transfer matrix SX associ- 
ated to the contact region in order to gain more knowl- 
edge and obtain more transparent expressions. We found 
the expression for the matrix current that describes a 
weak ferromagnetic contact to first order in SX (Eq.(|l9|)). 
In order to describe spin-flip processes at the contacts we 
need to go to higher orders in the asymmetry SX. For 
the case of a tunnel barrier, the tu nneli ng current takes 
a very simple and clear form (Eq.( 21a )). This expres- 
sion is characterized by three conductance parameters 
Gt ,Gmr and G$, which can be expressed in terms of 
the spin conductances G^q) and the (complex) mixing 
conductance G^. 

We thank discussions with Ya. M. Blanter, A. Brataas, 
and Gerrit E. W. Bauer. This work was financially sup- 
ported by the Stichting voor Fundamcntcel Onderzoek 
der Materie (FOM). 



APPENDIX A: PERTURBATION EXPANSION 



We perform a perturbation expansion of 
I + G\M< G 2 M) in terms of the parameter SX as 



fl + GiM * G 2 M) 1 = - 1 



KSM^- ] 



1 + .4 1 + A 1 + A 

(Ala) 



where 



h = Ml M 



(Alb) 



The circuit theory of mesoscopic transport is a system- 
atic way to describe transport in multiterminal hybrid 
structures in which general rules analog to the Kirch- 
hoff 's rules of classical circuits are used to solve the cir- 
cuit and compute the transport properties of the system 
under study. One of these rules imposes that the sum 
of all "matrix current" into a node must be zero. That 
is why, the expression of the "spin-charge" matrix cur- 
rent that flows between two nodes (or between a reservoir 
and a node) through a contact/ connector is as essential 
for the the circuit theory, as the Ohm's law is for clas- 
sical electric circuits. In this paper we have generalized 
the expression of such matrix currents for the case of 
multiterminal systems that included ferromagnetic and 
superconducting reservoirs connected through magneti- 
cally active contacts to one or several normal nodes. 

We have derived the most general expression for this 
matrix current Eq. ( |l4| ) in terms of of isotropic quasiclas- 
sical Green's functions Gi( 2 ) at both sides of the contact 



A = Q GiG 2 (Ale) 

SM^ = G^QqSX + GiSXiQo&a (Aid) 

To first order in SX an expansion in terms of normal 
eigenmodes is possible. The matrix M = (M) S ' B , 

can be diagonalized in the basis of eigenmodes N(N). 
For each eigenmode N, the elements of the transfer ma- 
trix (m)* , can be expressed in terms of spin-dependent 
transmission and reflection amplitudes at each side of the 
junction as 

m' +<+ = tj -2^ r i (V J r 2 (A2a) 

aB 

»>-r ( A2b ) 



G 



= t 



(A2c) 



(A2d) 



,where s,s/(a,/3) = {f, J.} are spin labels and 1(2) de- 
notes left (right) side of the contact region. Usually the 
transmission and reflection amplitudes are introduced in 
terms of the scattering matrix S 



S : 



r S S 



f-S s „s s 



-IS s 



(A3) 



where in this case t s s and r s s are matrices in Spin 
space. The scattering matrix S and the transfer ma- 
trix M are equivalent descriptions of a contact region. 
Nevertheless the transfer matrix M obeys a multiplica- 
tive composite rule, whereas the scatt er i ng matrix obeys 
a more complicated composition rule.EaEa. 

By writing in analogous way the other term in Eq. Jl4| ) 
(l — GiM' G*2 M) and after some algebra, we obtain the 
following expression for the matrix current / to first order 
in SX 



(A4a) 



where 



N,<j 



1 



i + g GiG 2 



(1 + 2 E 2 Gi -g GiG 2 ) 
(A4b) 



iW( S ) = -2e 2 Tr N , a U 



-0- 



with 



l + Q GiG 2 I + Q0G1G2 

i + e 2 Gi)} 



6 = [GtGiQaSX + Gi5X^Q G 2 



(A4c) 



(A4d) 



A general property of the transfer matrix is the flux 
conservation which can be expressed as 

E 2 M = E 2 (A5) 

Substituting Eq. @ in Eq.flA5|) we get 

(1 + SX*) M f E 2 Mq (1 + SX) = E 2 (A6) 

The matrix Mq should also obey flux conservation 



Ml E 2 M = E 2 



(A7) 



so we found the following equation for SX 

SX^ Z + E z SX + SX^ Z SX = 0. (A8) 



To first order in SX we have 

SX^ = — E 2 SX E 2 



(A9) 



From Eq.(AS) we obtain for the elements of SX 

SXl , = -SX+ + (AlOa) 



SX\ 



5X1 



SX^ 



-SX- 



(AlOb) 



(AlOc) 



We see that SX + + and SX are pure complex ele- 
ments and SX^ and SX |_ are complex conjugate. 

The matrix Mq 



Mo 



m 





0+ + 



5 



(All) 



is symmetric under time-reversal transformation, which 
involves both transformation in Spin and Channels space, 
under the following prescription 



a y t x (M )* £* a y = Mo 



(A12) 



being the x-Pauli matrix in directional space and a y 
the y-Pauli matrix in Spin space. As a result, the com- 
ponents of Mq fulfill the following relations 



,° t 



t 



(Al3a) 



(A13b) 



In general SX can be split in symmetric and anti- 
symmetric parts with respect to time-reversal 



'SX) =3 V H X (SX)*Y, x a y 



(SX) 



SX) = o y Y? (SX) t x a y = (SX) 



(Al4a) 



(A14b) 



(Al4c) 



Time-reversal transformation over Spin gives always an 
anti-symmetric contribution SX — S y (SX)* a y = — SX, 
so (fiX^j will be symmetric (anti-symmetric) with re- 
spect to transformation in directional space ( SX ) = 

\ ) A(S) 

t x (8X)*t x = ±(SX) A{S) . Taking into account the 
structure of SX given by the condition SX^ = —T, Z SX 
£ 2 , SX 



SX 



-SX_ _ SX+ _ 



SX* 



-SX: 



(Al5a) 



7 



we find for (sx) = (SX)* 2> = + (SX) 



A' 



SX = 



-5X_ 
5X1 



5X + 
-SX, 



SX^ 
5X1 



SX^ 
8X_ 



(A15b) 



In this case 5X + 
Space 



-SX and SXa is in directional 



sx A 



5X A 

sx: 



sx + _ 
-sx,. 



E*,E X , E v . (A15c) 



On the other hand for 8X S , (SXJ = Y, x (SX) S 1 



SX = 



-SX- _ sx + _ 
5X1 -5X+ 4. 



-5X1 



-5X+ _ 



(A15d) 



Now 5X + + = —SX , and 5X + _ must to be zero. As 

a result 5X$ is in directional Space 



SXc 



SX + + 
5X+ 



1. 



(A15e) 



For the matrix 5M — MqSX we have 
5M 



Sm^ 
5m. 



Sm^ 
8m,- 



(A16) 



T,a m + o 5X <y + J2a m + " 5X ° - 



where a = {+, — }. Analogy to M, the matrix S can be 
written 



S = S + 5S 



(A17) 



where Sq corresponds to the spin-independent part and 
5S accounts for the spin structure of S. Now we assume 
that 5S is a small deviation with respect to So (5S <C So). 
In this case the elements of S are 



= ri+ Sri 



t\ s =h + 5t 1 



r 2 s = r 2 + Sr 2 



(S s 

l 2 



h + 8t 2 . 



(Al8a) 



(A18b) 



(A18c) 



(A18d) 



The quantities t\(2) ( r i(2)) ar e spin- independent trans- 
mission(reflection) amplitudes, whereas 5t s x ^ (Srl^J 



are the spin-dependent correction to the total transmis- 
sion (reflection) t\ s fr| s J amplitudes. The elements of 
S can be expressed in terms of the elements of M like 



/ sis \ — 1 s\s 

-(m J ml 



iS s 

'2 



m,_(m /) 



tl s =(m*' s _)- 1 . 



(Al9a) 



(A19b) 



(Al9c) 



(A19d) 



By substituting Eq.( [A16| ) and Eqs. ( |Al8aHAl8d| ) , using 
analog relations between the elements of So and Mq and 
by expanding (m^T,/ - ) -1 up to order 6m , we obtain 

5r s 2 s ' = r 2 5X^ { - 5XJ ( + {r 2 5XJ: i - SXJ! L ) r 2 

(A20a) 



5t s 2 s =t 2 { 5X° S , +r 2 8X1 L ) 



5r{ 



ti h SX* s _ 



6t\ 



h(r 2 5X^1 -8XJL 



(A20b) 



(A20c) 



(A20d) 



So we see that the elements of SX can be expressed in 
terms of the transmis sion a nd r eflectio n amplitudes St 
(Sr). Note that in Eq. ( |A20a| )-Eq. ( |A20d| ) we explicitly in- 
clude the spin indices s, sf of the elements SX and St (8r) 
to emphasize that both SX and St (Sr), contain informa- 
tion about the spin structure of the contact region. 

Now we consider the case of SXa which corresponds 
to a weak ferromagnetic conta ct (^( 8X) A ~ a 

To evaluate the trace of Eq. ( A4c[ ) we can re- write the 
expression in the following way 



e 2 Trjv, CT {-2 



T 



N 



(A21) 



(G 2 Gi 



<i + T N ({G u G 2 } -2) 
T 1 ) (Gi&QoSX + Gi«tQ G 2 
. (1 + t z Gi) 



A + T N ({Gi,G 2 }-2) 



} 



The central term can be separate in four different terms 



as 



8 



QoSX G 2 Gi + G 2 5XtQ oGl + GiG 2 SX Qg 1 
+Gi Qq W G 2 



+ I (1) ( £ ) = e 2 



^4 + Tjv({Gi,G 2 } 



2) 



[{Sf N ,G 2 },G 1 ] (A29) 



(A22a) 



Qo^XQo 1 + G 2 SX^G 2 + GiG 2 SX G 2 Gi 
+G X Qo 1 ^ t Qo Gi 



(A22b) 



Q 5X G 2 Gi + G 2 5X^Q Gi + GiG 2 5A Q 1 \ ^ ■ 
+G 1 Qo 1 SX^G 2 Gl - 

(A22c) 



Qo^AQg 1 + G 2 £XtG 2 + GiG 2 SX G 2 Gi 
+Gi Qo W Qo Gi 



(S 2 G x ) . 
(A22d) 



The trace o ver "d irection" indices ct(ct') gives for the 
first term Eq.(|A22aD 



2Re[B* SX 12 ] G 2 Gi - GiG 2 (2Re[B*6X 12 }) (A23) 
+G 2 2 Re[_B* 5X i2 ] Gi - G x 2Re[5*5Xi 2 ] G 2 . 

Note that Gi(2) and 5Xi 2 are still matrices in "Keldysh- 
Nambu-spin" space and do not commute with each other. 
B* can be written in terms of reflection and transmission 
amplitudes like 



B 



2m + 



Ti 



-r 2 . 



N 



(A24) 



The spin-dependent corrections to the transmission and 
reflection probabilities are defined like 



#Tl(2) = t 1{2) 5*1(2) + ^t\ (2) t 1(2) 



(A25) 



4 + T w ({Gi,G 2 }-2) 

For the second term (Eq. |A22b[ ), using the prop- 
erty that the trace is cyclic, it is easy to see that 
Tr CT {SX} =Ti> {5At} =Tr a {Q 8XQ^} 

=Tr CT {Qq 1 5A t Qo} = 0, so the second term gives no 
contribution to I\. 

The third and fourth terms include the factor E 2 Gi 
Hin the right hand side. The presence of the matrix 
£ 2 will change the struct ure of this term in Channels 
space with respect to Eqs.A22a and A22b. Using that 
[SX, G 2 ] = 0, G 2 being the Green's function of a (nor- 
mal/ferromagnetic) reservoir, finally I\ gives 



E 

JV 



4 + T JV ({G 1 ,G 2 } -2) 



[{Sf N ,G 2 } + (4~2T N ) SO N 
2 

X 4 + T Ar ({G 1 ,G 2 }-2) 



4 5S iV ,Gi] (A30) 



5X + +i at and SQn — 
\ N both p ure i magina ry quan- 



being 5^, N = T N 
2 [5X+ + + i Im (r 2 SX+ _) 

tities. Finally inverting Eq. (|A20a| )-Eq. ( |A20d|) we can 
express also the quantities 5Sjv and 50jv in terms 
of spin-independent transmission(reflcction) ampli- 
tudes ii(2) (r 1 ( 2 )),and its the spin-dependent corrections 

5*1(2) ( Sr l(2)) 



55 



N 



se 



t\st 2 



N 



'•[Sri 



5t\ t 2 



T 



N 



(A31) 



(A32) 



5-Ri( 2 ) =rj (2) 5ri (2 ) + 5r| (2) ri (2 ) 
For this quantities we find 
5Ti 5T 2 SRi 8R 2 



(A26) 



Ti T 2 Ti T 2 

so 



= 2Re[r 2 5X^ S _]. (A27) 



2Re[B*5X + -'] = - (2 Re[r 2 5X + ^']) = -=-^. 

(A28) 



At first order in 5X, the contribution given by Eq.(A2£) 
is zero for in the case of [SXj s . Finally The contribution 
to this first term for the I\ current in terms of T/v, STn 
and Gi( 2 ) is 
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